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We study the possibility of realizing a growth rate of matter density perturbations lower than 
that in General Relativity. Using the approach of the effective field theory of modified gravity 
encompassing theories beyond Horndeski, we derive the effective gravitational coupling GeU and the 
gravitational slip parameter rj for perturbations deep inside the Hubble radius. In Horndeski theories 
we derive a necessary condition for achieving weak gravity associated with tensor perturbations, 
but this is not a sufficient condition due to the presence of a scalar-matter interaction that always 
enhances GeS- Beyond the Horndeski domain it is possible to realize Gefr smaller than Newton’s 
gravitational constant G, while the scalar and tensor perturbations satisfy no-ghost and stability 
conditions. We present a concrete dark energy scenario with varying ct and numerically study the 
evolution of perturbations to confront the model with the observations of redshift-space distortions 
and weak lensing. 


I. INTRODUCTION 

The observations of redshift-space distortions (RSD) 
and weak lensing CH , combined with Cosmic Microwave 
Background (CMB) measurements d, offer the possi¬ 
bility of testing General Relativity (GR) on cosmological 
scales. In particular, the observational evidence of late¬ 
time cosmic acceleration Q may be related to some mod¬ 
ification of gravity at large distances. The dark energy 
equation of state wde < ~1, which is allowed from the 
joint analysis of GMB and supernovae la (SN la) data 
[lOj, can be realized in modifie d g ravitational theories 
without ghosts and instabilities |lll |. 

If we modify gravity from GR, an extra scalar degree 
of freedom usually emerges due to the breaking of gauge 
symmetries of GR [l^. This scalar field mediates an 
extra gravitational force with a matter sector. In f{R) 
gravity, for example, the effective gravitational coupling 
Geff between the gravitational scalar and matter is 4/3 
times as large as Newton’s gravitational constant G in the 
regime where the scalar mass M is much smaller than the 
physical momentum /c/a of interest Il3ll. 

The recent observations of RSD [lj-[lg and cluster 
counts have measured the lower growth rate of mat¬ 
ter density perturbations Sm than that predicted by the 
A-Gold-Dark-Matter (ACDM) model. In fact, the Planck 
GMB measurements [T3, [3 ™ tension with the RSD 

data and the Hubble expansion data from SN la. One 
possibility for reconciling this discrepancy is to incorpo¬ 
rate massive neutrinos [l9| , but this increases the tension 
between the CMB and the Hubble expansion measure¬ 
ments Sim. 

Another possibility for realizing a lower cosmic growth 
rate is interacting models of vacuum energy and dark 
matter [2^-[3 (see also Ref. [13] )• If there is an energy 
transfer from dark matter to dark energy, it is possible to 
reduce the tension between the CMB and RSD measure¬ 
ments [3- Most of these interacting models are based 
on a phenomenological approach, in that the equations 
of motion do not follow from a concrete Lagrangian. In 


this case, even if a lower growth rate consistent with ob¬ 
servations is realized, it is not generally clear whether 
theoretically consistent conditions such as the absence of 
ghosts are satisfied or not. 

There exists a modified gravitational scenario-dubbed 
the Dvali-Gabadadze-Porrati (DGP) braneworld model 
[3~ which possesses an explicit Lagrangian in the five¬ 
dimensional bulk space-time. In the branch where the 
late-time cosmic acceleration occurs, it is known that the 
effective gravitational coupling Geff is smaller than G on 
scales relevant to large-scale structures [3. However, 
ghosts are present in this accelerating branch. Thus, in 
the DGP model, the lower cosmic growth rate is related 
to the appearance of ghosts [3lj| - 

Now, a question arises. Are there some modified 
gravity models with concrete Lagrangians realizing weak 
gravity (Geff smaller than G) on cosmological scales, 
while avoiding the ghosts and instabilities associated 
with the propagation speeds of scalar and tensor per¬ 
turbations? In order to address this problem, we fo¬ 
cus on a very general class of scalar-tensor theories 
dubbed Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theo¬ 
ries [HI- This also accommodates Horndeski theories 
[S^l^the most general scalar-tensor theories with second- 
order equations of motion in a general space-time. 

The action of GLPV theories has been derived in such 
a way that the Horndeski action written in the frame¬ 
work of the Arnowitt-Deser-Misner (ADM) formalism 
(3 does not obey two additional conditions. While 
this can generate derivatives higher than second order 
in a generic space-time, there is no extra propagating de¬ 
gree of freedom on a flat Friedmann-Lemaitre-Robertson- 
Walker (FLRW) background according to Hamiltonian 
analysis in terms of cosmological perturbations (3^ . ■ 

This conclusion also holds for odd-mode perturbations 
on a spherically symmetric background [37j| . 

In GLPV theories the tensor propagation speed 
squared can deviate from 1 even during the radiation- 
and early matter-dominated epochs, while in Horndeski 
theories this is limited by the two extra conditions men- 
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tioned above [s^, • Even for a simple canonical scalar 

field (f) with a potential, the deviation of from 1 can 
give rise to an interesting observational signature such as 
a large difference between two gravitational potentials 'll 
and $ [ 13 . For constant models one has Geff < G 
in the superluminal regime (c^ > 1), but Geff needs to 
be very close to G due to the fact that the scalar propa¬ 
gation speed squared Cg becomes negative as is away 
from 1. In this paper we show that it is possible to realize 
weak gravity in varying models, while satisfying the 
no-ghost and stability conditions associated with scalar 
and tensor perturbations. 

In Sec. m we begin with a brief review of the effective 
field theory of modified gravity in which GLPV 

theories are encompassed as a specific case. In the pres¬ 
ence of a matter component, we present the background 
and linear perturbation equations of motion in the uni¬ 
tary gauge in Sec. IIIII 

In Sec.|TV]we derive the effective gravitational coupling 
Geff/G and the gravitational slip parameter 77 = —d'/d' 
by employing a subhorizon approximation for the per¬ 
turbations relevant to large-scale structures. In GLPV 
theories time derivatives of metric perturbations are left 
even under this approximation, so the usual quasistatic 
approximation is trustable only when these time deriva¬ 
tives are suppressed relative to other terms. In this 
case, we implement such time derivatives as corrections 
to leading-order terms. 

In Sec. El we discuss the possibility of realizing weak 
gravity by expressing GeS /G and ry in terms of quantities 
associated with the no-ghost and stability conditions of 
scalar and tensor perturbations. In Horndeski theories 
we derive a necessary condition for realizing Geff < G, 
which is related with quantities appearing in the second- 
order action of tensor perturbations. However, this is not 
a sufhcient condition due to an extra scalar interaction 
with matter which always enhances Geff. 

In GLPV theories the value of is not restricted to 
be close to 1 even in the early cosmological epoch. In 
Sec. IV Bl we propose a simple model with a time-varying 
Cj in which the realization of weak gravity is possible 
without ghosts and Laplacian instabilities. In Sec. IVTl 
we numerically solve the full perturbation equations of 
motion for the decreasing cf model with 0 < c? < 1 
and show that the growth rate of matter perturbations 
associated with RSD measurements can be lower than 
that predicted by the AGDM model. 


II. MODIFIED GRAVITATIONAL THEORIES 


The most general scalar-tensor theories with second- 
order equations of motion are known as Horndeski theo¬ 
ries [ 3 ^ . The four-dimensional action of Horndeski the¬ 
ories is given by S' = / d^x^f—g L with the Lagrangian 



L = G2(<^,V)+G3(</),V)n<^ 

+G4(</>, X)R- 2G4.x(</>, X) 

+G5(</>, + ^G5,x{cl^,X)[{a^f 

-3(n</.) , (2.1) 

where g is a determinant of the four-dimensional met¬ 
ric gfj,u, G 2 , 3 , 4,5 are functions in terms of a scalar field 
(j) and its kinetic energy X = R and G^^ are 

the four-dimensional Ricci scalar and the Einstein ten¬ 
sor respectively, and a semicolon represents a covariant 
derivative with □(() = {g^'' 

GLPV theories [s^ correspond to the generalization 
of Horndeski theories derived by reformulating the La¬ 
grangian (12.1|) in terms of the 3-(-l ADM decomposition 
of space-time [b^ with the foliation of constant-time hy¬ 
persurfaces Ej. The ADM formalism is based upon the 
line element ds^ = g^^dx^dx'^ = —N'^dt^ + hij{dx'^ + 
N^dt){dx^ + N^dt), where N is the lapse, V* is the shift, 
and hij is the three-dimensional spatial metric. 

The extrinsic curvature and the intrinsic curvature 
are defined, respectively, by and = 

where = (—V, 0,0,0) is a normal vector 
orthogonal to St and is the three-dimensional 

Ricci tensor on St. In the following we shall focus on 
a flat FLRW background described by the line element 
ds^ = —dt^ + a?{t)5ijdx^dxG where a{t) is the scale fac¬ 
tor. In the unitary gauge the scalar field (j) depends on 
the time t alone and hence X = where a dot 

represents a derivative with respect to t. For this gauge 
choice, the action of GLPV theories can be written as 


5 = J d^x^/^L{N,K,S,n,U;t) 

J d xy/ g Lm{g (2.2) 

where Lm is the matter Lagrangian, and 
L = A2{N,t)+A3{N,t)K 

+Ai{N, -S) + Bi{N, t)n 

+A3{N,t)K3 + B3{N,t)(u-]^K'R^ . ( 2 . 3 ) 

Here we have defined K = S = R = 

U = and K3 = + 

\ with arbit rary functions A2,3,4,5 and i?4_5 
that depend on N and t [^[^. 

The coefficients G2,G3,G4,G5 in Horndeski theories 
are related to A2, A3, A4, A3 and i?4,H5, as Q 

A2 = G2 — XF3^^ , 

A 3 = 2{-Xf/‘^F3,x - 2y^G4,0 , 

A4 = —G4 -l- 2 VG 4 ,j 5 f -l- VG5^0/2 , 

B 4 = G 4 + X{G5,^-F3,^)/2, 

A5 = -(-V)3/2G5.x/3, 

B5 = -V^F 5 , 


(2.4) 
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where a comma in a lower index represents a partial 
derivative with respect to a given scalar quantity, and 
F 3 and F 5 are auxiliary functions satisfying 03=^3 + 
2XF3^x and = ^ 5 /( 2 X) + F^^x- From these re¬ 

lations it follows that Horndeski theories obey the two 
conditions 

Ai = 2XBi^x-B^, A3 = -^XB3 ,x. (2.5) 

GLPV theories are described by the Lagrangian (lO) 
without imposing the two constraints ( 1 ^ . Even with¬ 
out these restrictions, the linear perturbation equations 
of motion on the flat FLRW background remain of second 
order without having an extra propagating scalar degree 
of freedom [H, HI, [sg . 

For the study of growth of large-scale structures, we 
have taken into account the matter Lagrangian Lm in 
Eq. (12.21) . We assume that, in the frame described by the 
metric 5 ^ 1 / (dubbed the Jordan frame), there is no explicit 
coupling between the scalar field (f) and matter fields 
The matter energy-momentum tensor following from L^n 
is given by 


rp^v _ 2 (^(v^ 9 Lm) ^2 0 ^ 

In what follows we describe the matter component as a 
barotropic perfect fluid, whose background energy den¬ 
sity and pressure are given, respectively, by p and P. 


III. LINEAR PERTURBATIONS 

The expansion of the action (12.21) up to second order in 
the perturbations on the flat FLRW background gives rise 
to the background and linear perturbation equations of 
motion. This was already derived in Refs. [3J,[^I3, so 
we simply quote the results in the following. We consider 
four scalar metric perturbations A, ijj, E and tensor 
perturbations 7 ^- described by the line element 

ds^ = — (1 -I- 2A)dt^ + 2ditj}dtdx'' 

-|-a^(t)[(l -I- 2C)(5ij -I- 2didjE + jij\dx^dx^ , (3.1) 

and choose the unitary gauge 


J^ = 0, E = 0, (3.2) 

where 6(j} is the perturbation of (f). 

The background values (represented by an overbar) of 
the extrinsic and intrinsic curvatures are given, respec¬ 
tively, by and TZ^i, = 0, where F[ = d/a is 

the Hubble parameter. Then, it follows that K = 3iL, 
S = 3F[^, and H = U = 0. We consider the pertur¬ 
bations of these geometric scalars, e.g., 6N = N — 1, 
6K = K — 3H. The perturbations of the matter energy- 
momentum tensors are denoted as 

ST° = -6p, 6T° = d,6q, 6T; = SP6/, (3.3) 


where Latin indices correspond to components in a three- 
dimensional space-adapted basis. 

Expanding the action (12.21) up to first order in scalar 
perturbations and varying the first-order action with re¬ 
spect to 6N and 5a, we obtain the background equations 

L + L^n-3HF=p, (3.4) 

L-F -3HF = -P, (3.5) 

respectively, where F = L x + The matter com¬ 

ponent obeys the continuity equation 

p + 3H{p + P) = Q. (3.6) 

The second-order action for tensor perturbations de¬ 
rived from Eq. (IQ) reads 

= J d'^xa^qt5^^5^^ , (3.7) 

where 



and 


£ = L,ti + -^L,u + 2^L,u ■ 
We require the two conditions 


<?t > 0, 

c? > 0, 


(3.8) 

(3.9) 

(3.10) 

(3.11) 


to avoid the tensor ghost and small-scale Laplacian in¬ 
stabilities, respectively. 

(s) 

In the presence of matter, the second-order action S 2 
for scalar perturbations in GLPV theories was given in 
Ref. [3^,|3^,|4^. GLPV theories satisfy conditions for the 
absence of spatial derivatives higher than second order 
[33,[51|. Varying with respect to 5N, and (/, it 
follows that 


(2L,Ar -H L^nn - 6HW -(- 12L,sH^) 6N 

+ (3(-?F\w- 




W - iiv + £)^ = Sp, 


W5N -4L^sC = -Sq, 

y)+4(V + £)^ + ^d ( 
-3(p + P)SN = 36P, 
respectively, where 

W = L,kn + 2HL^sn + 4.L^sH , 
HL^xu, 


B = L^xn — 


V = -3Jg. 


(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 
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The continuity equations = 0 and ST^i.fj_ = 0 

lead, respectively, to 

dp + 3H {5p + 5P) 

= + (3.18) 

Sq + mSq = -{p + P)SN-SP. (3.19) 

Substituting Eq. (13.171) into Eq. (13.141) and employing 
Eq. (I3.19L we obtain 


{L,s + HL^s)'^ + +{V + £)6N + £( = 0 ■ (3.20) 


On using Eqs. (13.121) and (I3.13L the second-order 
f s) 

scalar action S 2 can be expressed in terms of C; the 
matter perturbations, and their derivatives [4^. Pro¬ 
vided that the matter component does not correspond 
to a ghost mode, the scalar ghost is absent under the 
condition [s^. IdM [5ll| 


9s = 


SqtilQqtWs + 3>V^) 


> 0 , 


(3.21) 


where Wg = 2 L,n + L,atat — dHW -f 12P['^L^s- Under the 
condition (13.101) . Eq. (13.211) translates to IQqtWs + SW^ > 

0 . 

In GLPV theories the scalar propagation speed 

a rared Cg is affected by the presence of matter [^, |^, 
. For nonrelativistic matter characterized by P = 0 
and 6P = 0, the value of Cg in the small-scale limit reads 


’-s 

9s 

where 


M + HM -£- 


\sP 


W2 


(1 + 2aH) 




V 


ttH = 


W 

£ 9 

- 1 = - 1 


^,5 


V 

L,s 


(3.22) 

(3.23) 

(3.24) 


The parameter an characterizes the deviation from Horn- 
deski theories . To avoid the small-scale instability of 
scalar perturbations, we require that 


c" >0. 


(3.25) 


The four conditions (j3.10p . ()3.1ip . (I3.2ip . and ()3.25p need 
to be satisfied for theoretical consistency. 

We define the gauge-invariant gravitational potentials 


^ = SN + ip, ^ = C + Hi}, (3.26) 


and the gravitational slip parameter 

$ 



Then, Eq. (13.201) can be written as 


(3.27) 


4 / + $ = J^(C-$)-(c2-l)C-aH5A^. (3.28) 

tlqt 


This shows that the deviation of 77 from 1 is induced by 
the variation of gt, the deviation of from 1, and the 
deviation parameter an from Horndeski theories. 

We also introduce the effective gravitational potential 

$2 = ^(4/-$), (3.29) 

which is associated with the deviation of light rays in 
weak leasing and CMB observations [s^. 


IV. SUBHORIZON PERTURBATIONS 

To study the growth of structures during the matter- 
dominated epoch, we shall take into account nonrelativis¬ 
tic matter satisfying P = 0 and SP = 0 for the La- 
grangian L^- We also define the gauge-invariant matter 
perturbation 

5m = 5-3Hv, (4.1) 

where 6 = 8p/p and v = dq/p. Taking the time derivative 
of Eq. (13.181) in Fourier space and using Eq. (I3.19p . we 
obtain 


^2 

dm + 2Hdm + ^4- = -3B - 6HB , (4.2) 

where fc is a comoving wave number and B = C, -\- P[v. 
The gravitational potential 4' works as a source term for 
the growth of matter perturbations. 

To estimate the evolution of 47, we consider the pertur¬ 
bations deep inside the Hubble radius, i.e., k/a:s> H. In 
Fourier space we employ the subhorizon approximation 
under which the dominant contributions to the perturba¬ 
tion equations are the terms involving dp and the terms 
multiplied by [6^. In GR, the accuracy of this 

approximation was numerically confirmed for the sub¬ 
horizon perturbations . 

By employing the subhorizon approximation in mod¬ 
ified gravity theories, we consider theories in which the 
deviation from GR is not so significant in a way that 
the terms other than those containing and dp are 

still subdominant to the perturbation equations, as in 
the case of GR. For example, the orders of the terms 
L,NdN and L,sH‘^dN in Eq. (13.121) are regarded to be at 
most of the orders of M^^H^dN. Under the subhorizon 
approximation, Eq. (13.121) reads 

^ [awX + 4(l-f Q!h)C] - (4.3) 



where 




W 

HL^s ’ 


X = H'lp, 


By] 


^ P 


(4.4) 


The perturbation y is related to the gravitational poten¬ 
tial 4>, as 4' = C-t-y. The definition of £tm comes from the 
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fact that the Friedmann equation (13.41) can be written as 

= p-A2- 6H^A5 - A2,n - ^HAs^n 
+6 H^L.sn + 12H^A5,n , (4.5) 

where we used the relation A 4 = —L g — 3HA^. Alterna¬ 
tively, one may introduce the General-Relativistic form 
of the matter density parameter tlm = p/{3H'^PI^{) by 

expressing Eq. (14.5|) in the form = p-\ -, where 

Mp\ is the reduced Planck mass. In this case, the rhs of 
Eq. (14.31) is replaced by L^g- 

Expressing Eq. (13.201) in terms of 4', (, x, and x, it 
follows that 

(1 -I- aH)4' + CtC + (1 + Cs - Q!Heif)x = , (4.6) 

where 


£5 = 


HL^g ’ 



(4.7) 


Under the subhorizon approximation the continuity 
equation (j3.18j) reads 

5p -H 3H5p ~ — iy-pi) + 4L,5C - yVSN^ , (4.8) 

where we employed Eq. (13.131) and ignored the term 
3C relative to {k'^/a^)'tjj■ Taking the time derivative of 
Eq. (j4.3p . we can eliminate the terms Sp and Sp on the 
Ihs of Eq. (14.81) . This process leads to 

CTvvll^ T (o^w T £w£^w T SHtti) X 

+4 [(1 + aH)(l + £ 5 ) + £aH] C — , (4.9) 

where 


W 

“ HW ’ 


= ^ 
H ■ 


(4.10) 


In the unitary gauge {Scj) = 0), the mass M of the 
scalar degree of freedom does not explicitly appear in the 
perturbation equations of motion. The above subhorizon 
approximation is valid for M smaller than Cgfc/a 
In the regime M ^ Csk/a, the scalar field is nearly frozen 
to recover the General-Relativistic behavior, so that GeS 
is very close to G. In fact, this happens for viable f{R) 
dark energy models in the early stage of the matter era 


Since an = 0 in Horndeski theories, the terms on the 
rhs of Eqs. (HU) and (HU) vanish. In this case, we can 
express x a-nd C as a function of d' from Eqs. (HU) and 
(14.91) and then derive 41 in terms of <5 by using Eq. (14.31) . 


In GLPV theories we need to deal with Eqs. (1461) and 
(HU) as the differential equations involving x and We 
define the quantities 


HC '' 


C-y - 


- X 
Hx 


(4.11) 


If and are smaller than the order of 1, the terms 
on the rhs of Eqs. (1461) and (HU can be regarded as 
the corrections to those on the Ihs. Let us consider this 
situation and express C and x in terms of and e^. 
Of course the evolution of and is known only by 
solving the full perturbation equations for a given theory, 
so Eqs. (HU, (HU, and HU are not closed for q;h ^ 0. 

We define the effective gravitational coupling GeS, as 


— ^ = -AnGeffpSm , (4-12) 


which works as a source term for the growth of on the 
Ihs of Eq. (14.21) . The potential (13.291) obeys 


^$s = -47rGEp(5^, (4.13) 

a‘‘ 


where 


1 -|- ry Geff 


(4.14) 


Recall that the gravitational slip parameter is given by 

V = -iC + x)/'^- 

It is convenient to express Gefi, rj, and E in terms of 
the quantities qt, Cj, Qs, and because the signs of them 
need to be positive. Substituting the relations £ = , 

V + £ = 4gt(l + an), W = AHqtaw, and M = 169t(l + 
a}i)/{Haw) into Eq. (I3.22|) . it follows that 


32(7t(l + an) 


l + 2 es-ey, + 

1 + an 

Cj. 6 H^ 1 -)- 2 crjr 

4(1 -I- an) aw 1 + an 


.(4.15) 


This relation can be employed to express Cw in terms of 
Cg. On using Eqs. (14.31) . (I4.6|) . and (14.91) with Eq. (14.111) 
and the relation 5m — 5 under the subhorizon approxi¬ 
mation, we obtain 


G 8 qt /2 ’ 


fl + /3 Ki 

2/2 8 qt ’ 


(4.16) 

(4.17) 

(4.18) 


where G = (87rAfpj) ^ is the gravitational constant, and 

























6 


/i = 64[2(1 + q;h)(1 - auf-H - + es)(l + + an + a^es + anec + Coh) + ^ auc^flmjqt 

—32cjavv [2 + e^n + + 2aH(l + Cs)] 9t + (^sCg + S^tC^), (4-19) 

/2 = (1 + Q^h){(1 + Q^h)<ZsCsCI^ + 8q;h[24(1 + — 4(1 + Q!h)q^w( 1 + £// + ~ ^c)]'?t} ’ (4-20) 

/a = (1 + aH){64[2(l + aH)(l + Cs + ctn + cknCs + ciHec 4" ^qh) + 3aHi2m]gt 

—32q!vv[1 + Cj + e^H + Es + Q;h( 2 + Cff + 2es + Ex)]*?* + ck^(qsCs + SgtCj)} . (4-21) 


Equation (|4.16l) shows that Geff is scale independent. 
This comes from the fact that we have ignored the field 
mass M relative to cjzja for its derivation. In other 
words, the results (I4.16p - (I4.18|) are valid in the regime 
M <C Csk/a. For the scalar degree of freedom associated 
with dark energy, the field mass is usually smaller than 
Csk/a at the late cosmological epoch. For some dark en¬ 
ergy models in which the chameleon mechanism is at 
work in the region of high density, there is a transition 
from the General-Relativistic regime (M ^ Csk/a) with 
Geff ~ G to the scalar-tensor regime (M <C Csk/a) with 
Geff given by Eq. (14.161) [H, . For the model discussed 

later in Sec. ED the condition M <C Csk/a is satisfied for 
subhorizon perturbations. 

In Horndeski theories (an = 0, = 0), the above an¬ 

alytic solutions are closed. In GLPV theories the terms 
auE^ and ohEx do not vanish in Eqs. (j4.19l) - (j4.2II) . In 
this case, we can check the validity of the subhorizon ap¬ 
proximation by solving the full perturbation equations of 
motion numerically for a given theory and by comparing 
the full results with the estimations (j4.16l) - (j4.18l) derived 
after the substitution of numerical values of and 
into Eqs. (14.19^ - 114.211) . In Sec. ED we shall do so for a 
concrete theory in the framework of GLPV theories. 


V. POSSIBILITY OF REALIZING WEAK 
GRAVITY IN THE COSMIC GROWTH HISTORY 

In this section we discuss the possibility of realiz¬ 
ing a gravitational interaction weaker than that in GR 
(Geff < G) on scales relevant to large-scale structures. 
We shall focus on GLPV theories described by the La- 
grangian (lOl) . Then, the tensor propagation speed 
squared is given by 


2 _ ^4 + B ^/2 

“ -A4- 3HA5 ■ 


(5.1) 


For the evaluation of the scalar propagation speed 
squared (I4.15p . it is convenient to express by using 
the background equations of motion (13.4p and ()3.5p . i.e., 


6 flm _ L.at + 

~ HW 


(5.2) 

I 


In what follows we shall discuss the cases of Horndeski 
and GLPV theories separately. 

A. Horndeski theories 


Substituting oh = 0 and e^jj = 0 into Eqs. (14.191) - 
(I4.21D , we obtain fi = a^gsC^c"^ -f 8(awCt - 4 - 4es)2gt, 
/2 = al^qsCs, and / a = + 8(q:w - 4)(awCt - 4 - 

4es)(7f Then, Eqs. (I4.16p - (l4.18p reduce to 


Geff ^ My, 

G 8qt 

_ alyqsC _ 

^ -k 8Q‘^qt 

8Q{Q -k Ow - 4) qt 
. + a?vc 2 (l + c?) 9s. 


pl(l+Ct^ 

Ibgt 


\ ai^cic/qsj 
2 -k 8(aw - 4:)Qqt 


(5.3) 

(5.4) 
,(5.5) 


where 


Q = ay - 4 - 4e5 . (5.6) 


Provided the rhs of Eq. (lOP is smaller than 1, it fol¬ 
lows that Geff < G. The contribution M‘^,c^/{8qt) in 
Geff/G originates from the tensor part, whereas the sec¬ 
ond term in the bracket of Eq. (15.31) comes from the in¬ 
teraction between the scalar field and matter. Under 
the no-ghost and stability requirements (13.101) . (I3.11L 
(13.211) . and (I3.25p . the latter contribution is always posi¬ 
tive. Hence the necessary condition for realizing a gravi¬ 
tational interaction weaker than that in GR reads 


My 

8 qt 


< 1 . 


(5.7) 


Due to the presence of the scalar-matter interaction, the 
condition (15.71) is not sufficient for realizing Geff < G. 
The quantities qt and are given, respectively, by 


q, — [A 4 + SHA^) , 

2 _ B4 + B5/2 

~ B4- 2XB4,x + HXB^^x ’ 


(5.8) 

(5.9) 


where we used the two relations (ESp. The explicit form 
of the quantity Q in Eq. (j5.6p is given by 


^3,AfCt + 477^.4 + 4A4 -k 477^.4(1 — C^) -k 677^^5 TvCj -k 12HA ^ -k 12[77^(1 — c /) -k 77)^5 

77(-A4 - 377 A 5 ) 


(5.10) 
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In GR we have —A 4 = B 4 = M^^/2, A 3 = 0, and A 3 = 
B 3 = 0, in which case qt = Mpj/ 8 , = 1 and Q = 
0. In modified gravitational theories, we generally have 
^piCtV( 8 gt) ^ 1 and Q 0 . 

As an example, let us consider theories described by 
the Lagrangian 


L = A 2 {N, t) + A 3 (A^, t)K + A 4 {t) (K^-S-n) , 

(5.11) 

in which case qt = —A 4 / 4 , = 1, and Q = —{A^^n + 

AA 4 )/{HA 4 ). The theory given by L = —e((/))X/2 — 
R(^) + M'^yF{(j))R/2, where e, V, and F are arbitrary 
functions of </>, belongs to a subclass of Eq. (15.1111 . For 
the choices e((/>) = ujEoMpi/c^ and F(r/)) = (j)/M„], this 
Lagrangian encompasses Brans-Dicke theory as a 
special case. On using the correspondence (12.4L the La¬ 
grangian of this theory is equivalent to Eq. (15.111) with 
the functions 

A3 = , 

A4 = -B4 = -^M^iF{cI,). (5.12) 


In this case the relation W = —T = Mpj(2i7F + (j)F^^) 
holds, so Eq. (14.151) reduces to Cg = 1 by using Eq. (15.211 . 
From Eqs. (I01)-(I53]1 it follows that 


Geff 


G 

F 





2Fe + 



Fe + 2Ml,Fy 


(5.13) 

(5.14) 

(5.15) 


The conditions qt > 0 and q^ > 0 translate to F > 0 
and 2Fe -I- > 0, respectively, so GeS is larger 

than G/F due to the presence of the second term in the 
bracket of Eq. (15.1311 . From Eq. (15.1411 the parameter 77 is 
smaller than 1. The enhancement of Geff is compensated 
by the smallness of 77 , so that we obtain the value (15.151) . 
In other words, we have Q + a-^ — 4 = 0 in Eq. (15.51) for 
the model discussed above. 

Let us consider the full Lagrangian (lOl) with the 
Horndeski relations ( 1 ^ . The behavior of the quanti¬ 
ties qt and Cj is crucial for the realization of the condi¬ 
tion (EZl). The quantity qt = L gjA is associated with 
the matter density parameter defined in Eq. (14.4L 
as flm = p/(24iL^(7t)- The matter perturbation equation 
(lO) can be expressed by using 17^, as 





qt \ 

“wCsC? 9s/ 


— 0 ■ 


(5.16) 

where a prime represents a derivative with respect to 
AC = In a. Note that we neglected the contribution on 
the rhs of Eq. (14.21) under the subhorizon approximation. 


In Eq. (15.161) the quantity qt appearing in the denom¬ 
inator of Eq. (15.31) has been absorbed into the definition 
of ilm. If ftm is smaller than that in GR due to large 
values of qt, it is possible to realize a cosmic growth 
rate smaller than that in GR. This is one possibility 
for the realization of weak gravity recently studied in 
Ref. [^. For this purpose we require that the scalar- 
matter coupling Q is suppressed to satisfy the condition 
Ct [ 1 -l-8Q^(7t/(Qf^CgCtgs)]llm < L In dark energy models 
based on f{R) theories, for example, we have that cf = 1, 
8 Q'^qt/{al^CgCtqs) = 1/3, and that the deviation of llm 
from the value of GR is not so significant, in which case 
the growth rate of 6 m is larger than that in GR [ilii. 

Let us proceed to the discussion of appearing in 
Eq. (15.161) . Due to the relations (12.5L the X depen¬ 
dence in B 4 and B 3 gives rise to the functions A 4 and 
A 5 depending on X as well as 4>. These terms contribute 
to the background equations of motion as the field en¬ 
ergy density and pressure. To realize a viable cosmol¬ 
ogy, the derivative terms 2 XB 4 ^x, HXB^^x, and 7 / 5/2 
in Eq. (15.91) should be suppressed relative to B 4 in the 
deep matter era, in which case is close to 1 . 

If the field (j) is responsible for dark energy, the deriva¬ 
tive terms mentioned above can be comparable to B 4 
after the onset of cosmic acceleration. This can lead to 
a deviation of from 1. Whether decreases or not 
depends on the models and initial conditions. In covari¬ 
ant Galileons [^, for example, the entry into the region 
Ct < 1 can occur for late-time tracking solutions [t^. 

From Eq. (15.101) the t and N dependence in A 4 and 
A 5 as well as the N dependence in A 3 lead to nonzero 
values of Q, which enhances the rhs of Eq. (|5.3D . Even if 
Cj starts to decrease after the end of the matter era, the 
scalar-matter interaction induced by the term Q should 
not be large enough to violate the condition Geff < G for 
the realization of weak gravity. In covariant Galileons, 
for example, this scalar-matter interaction usually gives 
rise to a Geff larger than G even for < 1 at the level 
of linear perturbations 0 . 

To summarize, the modification to Geff coming from 
tensor perturbations is the first crucial factor for realizing 
weak gravity in Horndeski theories, but the condition 
(EH) is not sufficient due to the presence of the scalar- 
matter interaction. In Sec. IV Bl we shall study what kind 
of difference arises in GLPV theories. 


B. Theories beyond Horndeski 

Since GLPV theories are not subject to the constraints 
(Ell), the values of are not restricted to be close to 1 
even in the early matter era. In GLPV theories, the 
general expressions of Geff, 77 , and S are not as simple 
as those in Horndeski theories due to the extra terms 
q;h and . To simplify the analysis, we focus on the 
theories described by the Lagrangian 

L = --X- V(/)) - + -^F{m , (5.17) 

























where F((j)) is a function of (j) different from 1. Unlike 
Ref. [ 3 ^, we do not restrict the situation to the case in 
which F{(j)) is constant. In fact, we will show that the 
theories with a time-varying F^cj)) allow for the realiza¬ 
tion of weak gravity in the regime 0 < <1. 

Since 2? = 0 for the Lagrangian (j5.17|) . it follows that 

an = - 1, c? = F {^). (5.18) 


dm, corresponds to H/F[^ = — 30^, where 

fix = is the density parameter of the field 

kinetic energy. Then, the scalar sound speed squared 
(I4.15P reduces to 


2 2 

Cs = Ct 


~ c^) + 2eo 
6 flx 


(5.19) 


We also have qt = = 0, W = 2i2Mpj = —F, 

aw = 4, and Cw = H /= —en- The background 
equation L x F T = which follows from Eqs. (1X1) and 


where = 2ctCt/i2. On using Eq. (I5.19p . we can ex¬ 
press (7s = <^^/(2i2^) without using (p^. Then, the quan¬ 
tities /i, / 2 , and /a in Eqs. (I4.19|) - (I4.21I) reduce to 


fi 

/2 

h 


8MpV? 


[2CgeQ,jj 2(Cg Ct)(Ct l)(Ct +£aH)£^^ 3(C(, l)c^^flrn\ 


-16M2jc2(c2 - 1) [{cl F F (c* - l){eH F 


Cs - C? 


[2cl{eH - encl F F c?(c? - l)(2eff - 311™)] - 16M^ict{cl - l){e^ - e ^), 

{2Cg[l -|- (cj — 2 — ch F Can) F ch ] — Ct (<3t ~ l)(2ct — 2 — 2eH F 311™ + 2eQjj)} 
Fl6M^icl{cl - l)(c?e,^ - e^). 


8M2 


plH 


(5.20) 

(5.21) 


(5.22) 


In the deep matter-dominated epoch characterized by 
Ch — 3/2 and H™ ~ 1, we take the quasistatic limits 
—>■ 0 and —>■ 0 in Eqs. (I5.20I1 - (I5.22E The validity 

I 


Geff 

G 


77 ~ 1-1- 


1-c? 5{l-cl){cl-cl) 
cl +c2c?(3-3c? + 2e„j"““’ 

5(1 - Ct)(Cs - Ct) 

3c?(l + c2 - cl) 3cl{cl -cl - l)[3c^ - 


When Ct is constant, i.e., Cm = 0, these results match 
with those derived in Ref. for scaling solutions re¬ 
alized by the potential U((/>) = If cl > 1, 

then Cg becomes negative for Ct — 1 > 2 flx/{flm — 211 x)- 
Provided that 1 < Cj < 1 -I- 211x/(llm — 2Hx), we have 
Geff < G from Eq. (15.231) . Since lljf <C 11™ during the 
matter-dominated epoch, we require that cl is very close 
to 1 to avoid the Laplacian instability of scalar pertur¬ 
bations. Then, the deviation of Geff from G is restricted 
to be small. 

If cl is smaller than 1 without any variation, it follows 
that Geff > G. In this case, for cl away from 1, cl tends 
to be much larger than 1 during the matter era. Hence 
Geff is also restricted to be close to G. In the regime 0 < 
Cj <C 1 , the parameter 77 is approximately given by 77 ~ 
1-|-5(1 —Ct)/(3ct) under the condition Cg » 1. This large 
deviation of 77 from 1 is the distinguished observational 


of this approximation will be checked in Sec. ED Then, 
the effective gravitational coupling (14.161) and the gravi¬ 
tational slip parameter (14.171) reduce, respectively, to 


(5.23) 


25cl(cl-cl)(l-cl) 


3ct (2 Fcl- Can) - 3(1 -b c2)c^(e«H - 1) + Screen] 


..(5.24) 


signature of the model (15.171) with cl smaller than 1 [s^ . 

If cl varies in time, both Geff and 77 are subject to 
change. If jean I is much smaller than 1 , we can ignore 
the contribution of the terms Ca^ appearing in the de¬ 
nominators of Eqs. (15.231) and (15.241) . Let us consider the 
case in which Cj < 1 and Cg is larger than the order of 1. 
Then, Eqs. (15.231) and (I5.24|) approximately reduce to 


Geff 

G 



77 ~ 1 -b 


5(1-eg) 

3cg 


25 

9c: 


,4 ■ 


(5.25) 

(5.26) 


The decrease (increase) of cg leads to the negative (pos¬ 
itive) values of = 2ctCt/F[. This means that, even 
in the subluminal regime (cg < 1 ), it is possible to have 
Geff < G for decreasing cg. For < 0, the parameter 
77 also gets increased. 














































9 


A model with varying can be realized for the func¬ 
tion 

F{(j)) = cl , (5.27) 

where ca and /3 are constants. For dark energy models in 
which V,x grows in time, Eq. (j5.19j) shows that [Cg | tends 
to increase as we go back to the past. This behavior can 
be avoided for sca ling dark energy models described by 
the field potential 

, (5.28) 

where Vd, V 2 , Ai, A 2 are constants with Ai > 10 and A 2 < 
1. During the matter era, the solution is in the scaling 
regime characterized by ilx = 3/(2Af) and Dm = 1 — 
S/Aj [t^. After the dominance of the second potential 
on the rhs of Eq. (15.2811 . the solution finally approaches 
an attractor with cosmic acceleration . 

We recall that Eqs. (I5.23|) and (|5.24l) are valid only 
during the deep matter era with negligible variations of 
C and X- In order to know the precise evolution of per¬ 
turbations from the matter era to today (especially after 
the onset of cosmic acceleration), we need to resort to 
numerical simulations. 



1 + Z 


Figure 1. Evolution of w^, cl, cl versus 1+z = 1/a for /3 = 0.1, 
cfi = 0.9, Ai = 10, A 2 = 0.5, and V^/Vi = 10“® with the 
initial conditions xi = ■\/6/(2Ai), X 2 = ■\/6/(2Ai), and 13 = 0 
at 2 = 735.4. The present epoch (the redshift 2 = 0) is 
identified by the condition flm = 0.3. 


VI. OBSERVABLES EOR A CONCRETE WEAK 
GRAVITY MODEL 


For the potential (I5.28P the scaling matter era corre¬ 
sponds to 


In this section we numerically integrate the linear per¬ 
turbation equations together with the background equa¬ 
tions for the theory given by the Lagrangian (15.171) . We 
consider the case in which the functions F((f)) and V((/)) 
are given, respectively, by Eqs. (j5.27l) and (15.281) . 


A. Background equations and propagation speeds 

In the presence of nonrelativistic matter the back¬ 
ground quantities xi = </> / {'/6H Mpi), = 

\/V/{'JZHMpx), and 2:3 = (ji/Mpi obey the equations of 
motion 

A = {xl -xl-l) + ^Ax^ , (6.1) 

^'2 = ^2:2 {xl - a;2 -I- 1) - ^Aa;ia;2 , (6.2) 

a;3 = VQxi , (6-3) 

where A = —MpiE^/V. The field equation of state 
= [(j)^/2 —V{(/))]/[4^/2+ V{(/))] and the matter density 
parameter Dm = p/{3H^M^l) can be expressed, respec¬ 
tively, as 


2 2 

^1 ~ ^2 O 1 2 2 

W(p — 2 , 2 ’ ^2 * 

XI + X 2 


(6.4) 


Xl = X 2 = 


2Ai 


W{p — 0 , Dm — 1 


A?’ 


(6.5) 


whereas the late-time scalar-field fixed point is charac¬ 
terized by 



X2 




Dm =0. 

( 6 . 6 ) 


From the big bang nucleosynthesis bound on Dm, the 
slope Ai is constrained to be Ai > 9.4 [75|. To realize 
the late-time cosmic acceleration we require that Af < 
2, under which the scalar-field-dominated fixed point is 
stable [tJ. 

In Fig. [ 1 ] we plot the evolution of for Ai = 10, 
A 2 = 0.5, V 2 /V 1 = 10“® with the initial conditions 


Xl = X 2 = \/6/(2Ai) and X 3 = 0. The field equation 
of state starts to evolve from = 0 and then it hnally 
approaches Wtj, = —0.917. Even if Xi and X 2 are initially 
away from these values, the solutions soon approach the 
scaling fixed point (16.5|) because it is a temporal attrac¬ 


tor [ 73 . After the dominance of the second potential on 
the rhs of Eq. (15.281) . the solutions are attracted by the 
fixed point (l(i.6l) . 

Since the scalar field evolves along the potential with 
velocity (/) > 0, the tensor propagation speed squared 
cl = decreases from the initial value cl for 

/3 > 0. Provided that 0 < < 1, always stays in the 
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subluminal regime (see Fig. [T] for the case = 0.9 and 
/3 = 0.1). 

Let us consider the stability condition associated with 
the scalar propagation speed given by Eq. (15.191) . First 
of all, the quantity can be expressed as 

tan = -‘iV&ficlxi , (6.7) 

which is negative for /3 > 0 and (/) > 0. Since the de¬ 
viation of Cj from 1 gives rise to negative values of Cg 
for Cj > 1, we focus on the case in which is in the 
subluminal regime. 

During the scaling matter era characterized by 
Eq. 1)6.51) , the condition Cg > 0 gives the following bound 


the perturbation equations of motion can be written as 


c'= <5iV + , (6.11) 

X' =(^^-^^X-cU5N + 0 , (6.12) 

<5' = + 3{xl - 1)6N + IC\ctC + K.) ,(6.13) 

V:, = -6N+^V^, (6.14) 

xlSN=^n„^Sm-^}C\c^,C + x), (6.15) 

where fC = k/{aH) and H'/H = —(3/2)(l + x\ — x^)- 


Taking the J\f derivative of Eq. (16.141) and using other 
equations of motion, we obtain 


/3< 


4Ai 


3 + (A? - 3) 



( 6 . 8 ) 


+ azFn = 



(6.16) 


When Ai = 10 and cl = 0.9, for example, /3 < 0.34. 

Substituting the values (16.61) of the scalar-field- 
dominated fixed point into Eq. (I5.19L it follows that 
Cg = Cj (1 — 4 /I/A 2 ). If we demand the condition Cg > 0 
for the future attractor, we obtain the bound 

/3<Y- (6.9) 

In the numerical simulation of Fig. [T] the parameter 
A 2 = 0.5 is chosen, so the condition dSH) translates to 
/3 < 0.125. We have numerically confirmed that, for 
(3 < 0.125, cl remains positive from the past to the 
asymptotic future. If we do not impose the condition 
Cg > 0 in the future, then the bound (16.91) is slightly 
relaxed, e.g., /3 < 0.15 for A 2 = 0.5. 

Figure [1] illustrates the evolution of Cg and cl for 
P = 0.1, Ctj = 0.9, Ai = 10, and A 2 = 0.5. Since cl 
decreases in time, this leads to the growth of Cg during 
the scaling matter era. Around the end of the scaling 
regime the ratio Dm/Dx starts to evolve toward 0, so Cg 
starts to the decrease. Finally, the ratio cl/cl approaches 
the asymptotic value 1 — 4/3/A2 = 0.2. 

We have thus clarified the range of /3 in which the 
stability conditions of scalar and tensor perturbations are 
satisfied. Since qt = A7pj/8 and = (jP there are 
no ghosts in our model. 


B. Perturbation equations and initial conditions 

Let us proceed to the discussion of perturbation equa¬ 
tions and their initial conditions. Introducing the dimen¬ 
sionless velocity potential 


( 6 . 10 ) 


where 


ai = 


4xi 


4A(1 — Xljn — xl) — •\/6xi(2 — Dm — ‘2x1) 1 

(6.17) 


«2 = (l-c?)^/C 
3 

2x1 


Dyj 

2x^’ 

3xi{4xi + 2(2Dm - 3)xi + Dm(Dm - 1)} 


'\/6A(Dm T 4x2^)(Dm 4“ Xi — 1) 


(6.18) 


The general solution to Eq. (16.161) can be expressed in 
the following form 




(6.19) 


where Vm'^ is the homogenous solution satisfying the dif¬ 
ferential equation Vm'' +aiVm'’ 4-02 = 0, and 14n ^ 

is the special solution. 

When cl < 1 the first term on the rhs of Eq. (16.181) 
can be much larger than 1 for subhorizon perturbations 
(/C 3> 1). During the scaling matter era characterized by 
Eq. (16.5p . the homogenous solution obeys 






9(A!-3) - 

2A? 


m 


0, (6.20) 


where clg = (Al — 3)(l — cl)/3 and the quantity /C evolves 
as lC{Af) oc If cl is constant, we obtain the follow¬ 

ing solution in the limit x = 2ceff/C 3> I : 



4- C 2 sin 



where ci and C 2 are constants. Thus, exhibits the 
oscillation whose frequency is related to Ceff. The ampli¬ 
tude of Vm'^ decreases in proportion to a~^. 


Vm = Hv, 
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1 + z 


Figure 2. Evolution of the perturbations (^,5N,Vm,X for 
the same model parameters and background initial condi¬ 
tions as those in Fig. [T] The initial conditions associated 
with perturbations are chosen to be = 1.1779 x 10“®, 
= -5.6540x10“®, 5 = 2.1253xl0“^ Vm = -2.6174x10“®, 
and A3 = 25 at redshift 2 ; = 735.4. 


Provided that is a constant smaller than 1, the 
special solution to Eq. (16.161) is approximately given by 
Vm^ ~ 2 x'ixl\VLrr,.(c^ — 1)] for subhorizon perturbations. 
In Ref. [39j it was found that the perturbations Id™, X, 
and C stay nearly constant along the special solution dur¬ 
ing the scaling matter era characterized by Eq. ira . 
This means that, as long as the homogeneous solution 
Vm^ is initially suppressed relative to the special solu- 
tion Vm , the latter remains the dominant contribution 
to Vm ■ 

For nonconstant , there is an additional term on the 
rhs of Eq. (16.161) . If we neglect the time derivatives of 

(s') 

Vm for the estimation of the special solution, we obtain 




Ol2 

2 


cl-I 


3 xi 


^1 1euH ^ 


( 6 . 22 ) 

(6.23) 


where, in the second line, we have picked up the first 
term on the rhs of Eq. ()6.18l) . From Eq. (16.231) we find 
that the change of cl actually gives rise to a variation 

of Vm'^ of the order of Vm'^ ~ —eQjjt4n’V(ct ~ !)• This 
implies that an oscillating mode may arise even for initial 
conditions of Vm that are very close to the value (I6.23|) . 

We recall that the parameter (3 is bounded from above, 
see Eqs. (16.81) and (16.91) . In particular the slope A 2 needs 
to be smaller than the order of 1 in order to be con¬ 
sistent with observations of cosmic acceleration, so the 


condition (16.91) gives the bound /? < 0(0.1). During 
the scaling matter era the parameter e^H i® given by 
etiH = ~6cj/3/Ai, so it is smaller than the order of 0.1. 

{g\> (g\" _ __ 

Then the terms aiVm and Vm in Eq. (16.161) are not 
as large as 02 for cl not very close to 1. 

Numerically we solve the perturbation equations of 
motion for the initial conditions Vm = Vm\ V^ = 0, 
and x^ = Oj where Vm^ is given by Eq. (16.221) . For given 
Sm, the initial values of 6 N^ ^ and x are known accord¬ 
ingly from Eqs. (16.121) . (16.141) . and (16.151) . The initial 
condition of Sm is chosen such that its value today is 
equivalent to (T8(0) = 0.82, where (T8(0) is the rms am¬ 
plitude of overdensity at the comoving 8 h~^ Mpc scale 
(h is the normalized Hubble parameter Hq = 100 h km 
sec“^Mpc“^). 

In Fig. we plot the evolution of Xj Kn, and 5N 
for the same model parameters as those given in Fig. [T] 
Unlike the constant cl model [s^, the perturbations I4i, 
C, and X do not stay constant even during the scaling 
matter era. This variation is induced by the change of cl 
appearing in the special solution (I6.23p . 

In Fig. [2]the perturbation 5N exhibits damped oscilla¬ 
tions with a non-negligible initial amplitude. This comes 
from the fact that SN is related to the derivative V^, 
as SN = —Vl^ -b {H'/H)Vm- We recall that, even by 
choosing an initial value of Vm that is identical to Vm ^ in 
Eq. (|6.22l) , the contribution to Vm from the homogenous 
solution Vm'’ arises due to the variation of cl. Taking the 
N derivative of Eq. (16.211) , it follows that the amplitude 
of Vm'’ is a;/2 times as large as that of Vm'’ ■ Since we 
are considering the case x^ 1, the oscillating mode par¬ 
ticularly manifests itself in the perturbation SN through 
the derivative Vm'’ ■ 

For initial conditions of Vm away from the special solu- 
(s) 

tion Vm , the perturbations exhibit damped oscillations 
with larger amplitudes than those plotted in Fig. [51 This 
situation is similar to what happens for the constant cl 
model [ 13 . Since it is likely that such large oscillations 
can be severely constrained from CMB observations, we 
focus on the case in which Vm is initially close to Vm in 
the following discussion. 


C. Observables 

We study the evolution of observables associated with 
RSD, weak lensing, and CMB. In Fig. [3] we plot the 
gauge-invariant gravitational potentials —'F and 4) as well 
as — 4 )e versus 1 -I- z for the same model parameters and 
initial conditions as those used in Figs. [I] and |5J Since 
'F = SN + x' ~ {H'/H)x, the oscillation of SN seen in 
Fig. [pleads to that of —'F. In Fig.|3]the oscillating ampli¬ 
tude of —'F is not so large due to the choice of the initial 
condition \Vm'’\ "C |14n^|- With the decrease of Vm'’ , 
the oscillations of —'F damp away. The gravitational po¬ 
tential —'F exhibits an overall decrease from the matter 
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Figure 3. Evolution of the gravitational potentials —'I', 
— <1>E for the same model parameters and initial conditions as 
those given in the captions of Figs.[T]and[21 


era to today, so this case corresponds to weak gravity. 

Figure S] shows the evolution of Geff/G computed from 
the definition of Geff given in Eq. (14.121) . From the scaling 
matter era to today we have Ggff/G < 1, while satisfying 
the no-ghost and stability conditions of tensor and scalar 
perturbations (see Fig. [T}. The realization of weak grav¬ 
ity comes from the fact that the last term on the rhs 
of Eq. (15.231) is negative due to the decrease of . For 
constant Cj smaller than 1, GeS is larger than G for sub¬ 
horizon perturbations. 

In Fig. m we also plot Ggff/G given in Eq. (14.161) derived 
under the subhorizon approximation. We have computed 
the terms and by solving the full perturbation equa¬ 
tions and then substituted them into Eqs. (14.191) and 
(14.201) . The numerical simulation of Eig. 2] shows that 
the approximate formula (14.161) can reproduce the nu¬ 
merical results with high accuracy for the modes deep 
inside the Hubble radius. Even by setting = 0 = in 
Eqs. (14.191) and (|4.20l) we find that the analytic formula of 
Ggff/G is a good approximation during the scaling mat¬ 
ter era, but it starts to deviate from the full numerical 
integration around the end of the matter era. In GLPV 
theories the quasistatic approximation ignoring the vari¬ 
ation of and is not trustable at the late cosmological 
epoch. 

In the regime where is larger than the order of 1, the 
gravitational slip parameter is given by Eq. (15.261) under 
the subhorizon approximation in the deep matter era. 
Since the second and third terms on the rhs of Eq. (15.261) 
are positive for < 1 and < 0, r; is larger than 
1. In Fig. 2] we find that ry grows as the decrease of 
from the matter era to today and it starts to decrease 



Figure 4. Evolution of Geft/G and 77 (black solid lines) for the 
same model parameters and initial conditions as those given in 
the captions of Figs. [T]and[3 The dashed red and blue lines 
correspond to the evolution of Ges/G and ry, respectively, 
derived under the subhorizon approximation involving and 
tx, he-, Eqs. (14.161) and (14.171) . 


in the future. This deviation of ry from 1 is one of the 
distinguishing features of our model. As we see in Fig. 21 
the subhorizon approximation based on Eq. (I4.17|) with 
inclusion of and is a trustable prescription for the 
evolution of ry. 

Since ry is larger than 1, the two gravitational poten¬ 
tials $ and — T are different from each other. In Fig. [3] 
we Hnd that <I> grows during the matter era, while —'F 
decreases. We recall that the weak leasing gravitational 
potential <i)i; obeys Eq. (14.131) with S given by Eq. (14.141) . 
Since Geff < G and ry > 1 during the matter era in the 
numerical simulation of Eigs. [3]and21 the small value of 
Geff appearing in S is compensated by the large value of 
ry. Hence —4)i; does not vary much relative to $ and —T 
(see Fig. 2]). 

The growth rate of matter density perturbations 6 m 
can be measured by peculiar velocities of galaxies in RSD 
surveys [zSlzl- Usually, this is quantified by the data of 
f{z)as{z) at redshift z, where / = 6 m/{H6m)- In Fig. [5] 
we plot the ten data points of f{z)as{z) with error bars 
derived from the measurements of 2dFGRS [I| , 6dFGRS 
0, WiggleZ 0, SDSSLRG 0, BOSSCMASS @, and 
VIPERS 0. The latest Planck measurement of the GMB 
power spectra provided the bound crg(O) = 0.829 ± 0.014 
at the 68 % conhdence level [2l[. 

In Fig. [ 5 ] we also show the theoretical curves for 
/3 = 0,0.05,0.1,0.12, c?. = 0.9, and cr8(0) = 0.82 
with the initial conditions of perturbations satisfying 
\Vm \ <C |14n^|. When /3 = 0 we have = 0 and 
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Figure 5. Evolution of fag in the redshift range 0 < 2 < 1 for 
c?i = 0.9, Ai = 10, A 2 = 0.5, Vi/Fi = 10“®, (78(0) = 0.82, and 
/3 = 0,0.05, 0.1, 0.12 (from top to bottom) with the initial con¬ 
ditions xi = v^/(2Ai), X 2 = \/6/(2Ai), xg = 0, and /C = 25 
ai z = 735.4. The initial conditions of perturbations are cho¬ 
sen to satisfy \Vm'^\ <C The bold dotted curve corre¬ 

sponds to the evolution of fag in the ACDM model. The data 
points with error bars come from the recent observations of 
f{z)ag(z) by 6dFGRS, 2dFGRS, WiggleZ, SDSSLRG, BOSS- 
CMASS, and VIPERS surveys. 


Geff/G ~ 1 -I- (1 — Cj )/Cs in the deep matter era, but Ges 
is close to G due to the condition » 1 [s^. In Fig. [5] 
the theoretical curves for /3 = 0 and the ACDM model 
do not fit the recent RSD data very well. 

When /3 > 0 the effective gravitational coupling Geff 
is smaller than G from the matter era to today, so the 
growth rate of 5m gets smaller relative to the case /3 = 0. 
For larger (3, the theoretical curves shift toward smaller 
values of fag. From Fig. [5] we find that the varying 
models with /3 > 0.05 show a better agreement with the 
RSD data compared to the constant model. We recall 
that /3 is constrained as /3 < A 2/4 to avoid Cg < 0 for 
the future attractor, i.e., /3 < 0.125 in the numerical 
simulation of Fig. [51 In the regime 0.05 < (3 < 0.125, the 
models are compatible with the RSD measurements for 
the values of (T8(0) constrained by Planck. 


VII. CONCLUSIONS 

In this paper we have studied the possibility of realiz¬ 
ing a gravitational interaction weaker than that in GR on 
scales relevant to large-scale structures and weak lensing. 
We have employed the approach of the effective field the¬ 
ory of modified gravity encompassing both Horndeski and 
GLPV theories as specific cases. The important quanti¬ 


ties associated with the no-ghost and stability conditions 
of tensor and scalar perturbations are given by Eqs. (EH), 
(13.2111 . (13.2211 . All of the quantities qt, c^, qs, and Cg are 
required to be positive for theoretical consistency. 

Since our interest is the evolution of perturbations for 
modes deep inside the Hubble radius (fc/a H), we have 
exploited the subhorizon approximation under which the 
dominant contributions to the perturbation equations are 
those involving and the matter density perturba¬ 

tion 5p. In Sec. IIVI we have derived the general expres¬ 
sions of the effective gravitational coupling Geff, gravita¬ 
tional slip parameter 77 , and weak lensing parameter E 
by using the quantities qt, c^, qs, and Cg. In GLPV theo¬ 
ries the time derivatives C and x do not vanish even after 
employing the subhorizon approximation, so we need to 
know the numerical values of and in Eqs. (14.1911 - 
(14.2111 for the computations of Geff, rj, and E (unless |e^| 
and I | are much smaller than unity). 

In Horndeski theories the analytic expressions of Geff, 
77 , and E are of the simple forms (I5.3II - (I5.5|I . In this 
case the necessary condition for the realization of weak 
gravity is given by Eq. ([53, but this is not a sufficient 
condition due to the presence of additional scalar-matter 
interactions. The scalar-matter coupling, which always 
enhances Geff, should not be so large as to give rise to 
values of Geff larger than G. 

In GLPV theories the two conditions (I2.5|l are absent, 
so Cj can deviate from 1 even in the deep matter era. We 
have presented a simple time-varying model described 
by the Lagrangian (15.1711 with the function (j5.27p . In this 
model the scalar propagation speed squared Cg, which is 
given by Eq. (15.1911 , increases to be much larger than 1 for 
dark energy models with decreasing fix toward the past. 
For scaling models described by the field potential (15.2811 
the ratio flm/flx stays constant during the matter era, 
so the unbound growth of Cg in the past can be avoided. 

The analytic expression of Geff /G given in Eq. (15.2511 , 
which is valid in the deep matter era with Cg ^ 1 and 
|e(^| -C 1 , [fix I shows that it is possible to realize weak 
gravity (Geff < G) for the decreasing cl model (Can 0) 
in the regime 0 < Cj < 1. In this case, the parameter 77 
is larger than 1 from Eq. (|5.26ll . 

In Sec. EH we have numerically solved the full pertur¬ 
bation equations of motion for the model (j5.17p with the 
functions (15.2711 and (j5.28p . The solution to the normal¬ 
ized velocity potential Vm = Hv can be written by the 
sum of the homogenous solution Vm'^ and the special so¬ 
lution Vm'^. Provided that Vm'^ is initially much smaller 

is) 

than Vm , the amplitudes of oscillating modes of pertur¬ 
bations are suppressed apart from 5N and that are 
related to the time derivative of Vm'’■ The numerical 
simulations of Figs. [3] and [4] show the realization of weak 
gravity from the matter era to today with 77 larger than 
1 . 

In Fig. m we have also computed the evolution of fag 
in the redshift range 0 < z < 1 for several different val¬ 
ues of (3. For larger f3 the theoretical values of fag get 
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smaller, so that these cases exhibit better compatibility 
with the recent RSD data relative to the case /3 = 0. 
It remains to be seen whether or not the future RSD 
data combined with other observational probes favor the 
lower growth rate of matter perturbations than that in 
the ACDM model. 

There are several issues we have not addressed in this 
paper. First, under the so-called disformal transforma¬ 
tion [78l - [80l| . the Lagrangian (I5.17|) can be transformed 
to the one in the Einstein frame in which the tensor prop¬ 
agation speed is 1 [H, [HI HH ■ Since a nontrivial kinetic 
coupling with the scalar field and matter arises in the 
Einstein frame [H, |H, js^ , the role of such an interaction 
should be understood in view of a coupled dark energy 
and dark matter scenario [ssj) . 

Moreover, it will be of interest to study the screening 
mechanism of the fifth force [s^ in local regions of the 


Universe for the model (15.171) . In some modified gravity 
models like the quartic Galileon, it was shown that the 
screening mechanism can give rise to a GeS smaller than 
G in the nonlinear regime of matter perturbations [s^ • It 
will be of interest to see whether such properties persist 
in more general modified gravity models in the framework 
of GLPV theories. 
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